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Noise correlations, such as those observable in the time of flight images of a released cloud, are calculated 
for hard-core bosonic (HCB) atoms. We find that the standard mapping of HCB systems onto spin- 1/2 XY 
models fails in application to computation of noise correlations. This is due to the contribution of multiply 
occupied virtual states to noise correlations in bosonic systems. Such states do not exist in spin models. We use 
these correlations to explore quantum coherence of the ground states and re-address the relationship between 
the peaks present in noise correlation and the Mott phase. Our analysis points to distinctive new experimental 
signatures of the Mott phase. The importance of these correlations is illustrated in an example of a quasiperiodic 
potential that exhibits a localization transition. In this case, in contrast to the momentum distribution, the noise 
correlations reveal the presence of quasiperiodic order in the localized phase. 



In recent years, great experimental progress has been 
achieved in the coherent control of ultra-cold gases. In partic- 
ular by loading a Bose-Einstein condensate into a tight two- 
dimensional optical lattice, an array of one dimensional tubes 
has been created |3l |4j |5j |a CD. Using this setup, recent 
experiments have been able to successfully enter the Tonks- 
Girardeau regime (TG) 01 El 01- where the strong interac- 
tions between the bosons mimic the Pauli exclusion princi- 
ple 01, and to realize the Mott insulator transition in one 
dimension @|. 

Recent theoretical and experimental studies have shown 
that the atomic shot noise |9l |lfj, ll 211 in the time of flight 
images can be useful to decode various correlations underly- 
ing many-body states of trapped ultra-cold atoms. For ultra- 
cold bosons in a Mott insulator state, for example, these sec- 
ond order correlations have been proven to complement the 
standard, momentum-distribution based, characterization of 
the phase coherence. . Such correlations have also been 

used to probe pair condensation in a fermionic superfluid 1 10]. 

In this paper, we develop a theoretical framework to com- 
pute noise correlations in one dimensional hard core bosonic 
atoms (HCB) in a lattice. The HCB Hamiltonian is identical 
to that of spin- 1/2 XY model which in turn can be mapped to 
that of spinless fermions. ll3l 1 1411 Although these three sys- 
tems have identical spectra and local observable, their off- 
diagonal correlation functions differ. Experimentally rele- 
vant two-point correlation functions are identical for HCB and 
spin- 1/2 systems, and the general formulation to calculate 
these functions was developed by Lieb and Mattis| 14], based 
on Wick's theorem. To the best of our knowledge, this for- 
mulation has not been extended to higher order correlations. 
Here we complete such an extension to treat four-point corre- 
lation functions. One of the central results of this paper is the 
discovery of important differences between higher order cor- 
relation functions of HCB and spin- 1/2 systems. The root of 
this difference is the fact that HCB systems have virtual states 
which may be multiply occupied, whereas spin- 1/2 states are 
always at most singly occupied. 

We use the noise correlations to explore quantum coher- 
ence in systems with and without an external parabolic con- 
finement. Our analysis generalizes, to the strongly correlated 



(fermionized) regime, previous studies of noise correlations 
carried out in the Mott insulator limit.l9l ll ill We show that, 
independent of the filling factor of the system, HCB exhibits 
second order coherence displayed as peaks in the noise shot 
images which reflect the order induced by the lattice poten- 
tial. This suggests that the second order coherence in noise 
correlations is a generic attribute of the strongly correlated 
regime and does not merely indicate reduced number fluctu- 
ations (This was also noted in earlier studies in a different 
context. |0 [3)- On the other hand, we find that the intensity 
of peaks depends exclusively on the relative coordinates only 
in Mott-insulator limit. Thus a regular pattern in the noise cor- 
relation could serve as a definitive signature of the existence 
of a Mott phase. Finally, we show that noise correlations are 
also important in describing the quantum coherence of HCB 
systems in the presence of quasiperiodic disorder. Quasiperi- 
odicity induces an Anderson-type localization transition iTnll 
in the system. In the delocalized phase, both the first and the 
second order correlations are found to exhibit characteristic 
Bragg peaks at the Fibonacci sites reflecting quasiperiodic or- 
der. However, in the localized phase only noise correlations 
show these peaks. 

The Bose-Hubbard Hamiltonian describes bosons in optical 
lattices when the lattice is loaded in such a way that only the 
lowest vibrational level of each lattice site is occupied and 
tunneling occurs only between nearest-neighbor sites, 1 15]: 

H = - J J2 a\a 3 + h ^ - !) + Yl V i h i C 1 ) 

Here a,j is the bosonic annihilation operator of a particle at 
site j, hj — a'-ctj, and the sum (i, j) is over nearest neighbors. 
The hopping parameter J, and the on-site interaction energy 
U are functions of the lattice depth. Vj represents any other 
external potential such as a parabolic confinement or on-site 
disorder. 

In a typical experiment, atoms are released by turning off 
the external potentials at time t = 0. The atomic cloud ex- 
pands , and is photographed after it enters the ballistic regime. 
Assuming that the atoms are noninteracting from the time of 
release, properties of the initial state can be inferred from the 
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spatial images: (h[x(t)]), reflects the initial momentum distri- 
bution, n q , and the image shot noise, Q[x(t), x'(t)] reflects the 
momentum space fluctuations, namely the noise correlations, 
A(?i,g 2 ), 



(n[x(t)]) oc (n q ) =J2 e 



g[x(t),x'(t)] oc {(n qi - {n qi ))(n qi - (h qi ))) = A(qi,q 2 ) 

n,m,l,j 

where L is the number of lattice sites. In the strongly corre- 
lated regime, Eq. (1) can be replaced by the HCB Hamilto- 
nian, 

H J^(5jS j+1 + b] +1 bj) + ]T Vi*j (2) 



Here bj is the annihilation operator at the lattice site j which 
satisfies [b^j, Sj] = 0, and the on-site condition b 2 = b^ 2 = 
0, which suppresses multiple occupancy of lattice sites. The 
same relations are fulfilled by Spin- 1/2 raising and lowering 
operators. However, the exact on-site commutation relation 
differ between HCB and spin- 1/2 operators, and this becomes 
important in processes involving virtual states. This distinc- 
tion between spin and HCB models has not been explored in 
earlier studies; it has no effect on either the mapping of HCB 
to free fermions or on the calculations of local observable and 
the momentum distribution. However, the presence of multi- 
ply occupied lattice sites in virtual states that occurs in HCB 
strongly affects the noise correlation. A simple example of 
this can be seen in the computation of the correlation func- 
tion, (l|6&t|l). In the spin- 1/2 model this correlation func- 
tion is zero, while for HCB it is equal to 2. On the other hand 
(l|&t&|l) = 1 for both systems. 

We now describe the way in which we have generalized the 
approach of Lieb and Mattis. We have found a simple recipe 
to take into account the problem of multiple occupancy of the 
virtual state. In our calculation of the four point correlations 
{b^bmbjbj), each occurrence of a term bjbj is replaced by the 

1 + b\b,j (Note that no new rules are needed if we encounter 

number operators, b^bj or operators at different sites). In other 
words, only those four-point correlation functions involving 
two or more equal sites where the pair has the form bjbj, 
should be treated differently from that of the corresponding 
spin correlations. We will refer this recipe as multiple occu- 
pancy of virtual state rule (MOV). The validity of MOV was 
checked by comparing various correlation functions obtained 
using the above recipe with those obtained by diagonalizing a 
full Bose-Hubbard Hamiltonian with a large U value. 

The following procedure is used to calculate the four-point 
correlation functions. We first rearrange the operators so that 
the site index is ordered (this is only relevant for the case 
when three or more site indices are different); then we ap- 
ply MOV; next, we use the Jordan- Wigner transformation 



accordingly to the prescription of Lieb and Mattis; and fi- 
nally we use Wick's theorem to write higher order corre- 
lations in terms of the free-fermionic propagators, gi m = 
SfL"^ 1 V'* V'm > with N being the total number of atoms 
and ipf^ the s th eigenfunctions of the single-particle Hamil- 
tonian -J(^ ( +i + Vii/jI s) = E^i/j^. To present our 
results we denote the creation and the annihilation operators 
by b a , where a = +1(— 1) for annihilation (creation) opera- 
tors, respectively, and the "site ordered" four-point correlation 
function is designated by Xabcd ■ 



blbmbjbj) 



X-abcd 



(3) 



In this equation it is implicit that a < b < c < d, and 
this order is implied in all the expressions that follow. We de- 
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Sij and Bj 



= (b^bj). The latter can be 
We introduce four matrices 



calculated in terms of Gij 11411 . 
M, S, X and Y in terms of which our results for the correla- 
tion functions can be written. These matrices are presented at 
the end of this paper due to their notational complexity. Our 
results for the correlation functions then take the form: 



X-abbd 



Xabcc 



X-abcd 
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i|M(a,6,d)|+Q+«/j 1 -, ) B,„, 
~|S(a,c,o0| + Q+<5 Q ,-i ) LI.,, 
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X(a, 6, c, d)\- 
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(<5 7 ,-i - <5 Q ,-i)|Y(a, b, c, d)\ 



(4) 



When three, four or two pairs of indices are equal the cal- 
culation is rather straight forward. The non-vanishing corre- 
lation functions of this type are given by: (6^ b n b' n b n ) = g nn , 

(bnbnb m bm) — 9nn9mm ~ 9nm> {b nb\nb m b n) = g n n9mm ~ 

9n m + 9nn and (SjjS„&t b m ) = $\J)J)\fi n ) = B nm with 
n ^= m. 

The formulas above were used in numerical calculations 
of the momentum distributions and the noise correlations for 
lattices of L = 55 sites and N = 19 atoms with periodic 
boundary conditions. Fig.l shows both first and second or- 
der correlations along with the density profile . There we also 
include effects of a magnetic confinement Vj = j 2 Q/J on 
various correlations. Experimentally the ratio il/J can be 
changed either by changing the lattice depth or the external 
magnetic confinement. The parameters used in our analysis 
relate to typical experimental set-ups such as the ones reported 
in ref |31 • 

The density profiles (Fig la) for different values of Q/J 
show that only in the case fl/J — 0.17 does the ground state 
of the system correspond to a Mott insulator. In this case 
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all the central N sites have unit filling. For Q,/ J = 0.018 
localization takes place only at the central site, while for 
il/J = 0.008 and 0, all the sites have filling factor less than 
unity. The formation of a Mott state with reduced number 
fluctuations for fl/J — 0.17 is clearly signaled by the mo- 
mentum distribution (Fig lb), which shows a flat profile. On 
the other hand for the other cases where number fluctuations 
are important, there is a clear peak in the momentum distribu- 
tion at q = 0. This peak reflects the quasi-long-range correla- 
tions in the density matrix of these systems I I 1611 . 

In contrast to the momentum distribution, the noise- 
correlations A(qi , 0) (Figlc) for all values of Q,/ J are peaked 
at qi = 0. This lattice induced peak is observable due to the 
strong particle correlations in the fermionized regime and dis- 
appears in the weakly correlated regime where most of the 
atoms are Bose-condensed. In our numerical study without 
any confining potential, where we vary the filling factors (not 
shown in the figures), the central (qi — 0) peak in the noise 
correlations showed the same particle-hole symmetry as the 
central peak in the momentum distribution. The value of 
A(0,0) increases with filling factor, v, up to v = 1/2, and 
then decreases to its minimum value of 2 — 1/L at unit fill- 
ing when the system is in the Mott phase. For partially filled 
states (N < L), the quantity A(qi, 0) has a dip near q\ = 0. 
This satellite dip disappears in the Mott phase, in which case, 
A(gi ^ 0, 0) is a constant equal to —1/L. As seen in Fig. 1, 
these observations remain qualitatively valid in the presence 
of a trap, where the analysis of the noise correlations is more 
complicated, 
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FIG. 1: Density (a), momentum distribution (b) and noise correlation 
(c) with different trapping potentials: fl/J = O(black), 0.008(blue), 
0.018(green), 0.17(red)). Here q 2 = 0, N = 19 and L = 55. 
For finite fi, the correlation functions are renormalized by a scaling 
factor N/Z where Z are the number of sites with non-zero density. 




FIG. 2: Mott (a) and non-Mott(b) phase (obtained with fi/J = 0.17 
and fi/J = 0.008 respectively) using L = 55 and N = 19. Each 
color corresponds to A (91,92) for a fixed 92 as a function of qi. 
Thus different colors correspond to different values of 172. Peaks 
occur when qi = 92 • Red (blue) show the peak when 172 = in the 
Mott ( non-Mott) phase. Other peaks correspond to 172 = 15(black), 
30(green), 45(yellow). Panel (c) is a blowup of panel (b) where we 
show (72 7^ correlations only. 

The existence of second order coherence in HCB systems, 
independent of their filling factor, implies that the peaks can- 
not be used as a signature of the Mott insulator in ID systems. 
Nevertheless, our numerical calculations show that only when 
the system is a Mott insulator, the noise-correlations exhibit a 
regular pattern, i.e. A ((ft, 52) ~ A(gi —92) (small differences 
seen in Fig. 2 are due to the finite trap). Fig. 2 illustrates the 
contrast between the Mott (Fig2a) and the non-Mott (Fig2b-c) 
phase. In addition to the variation in the intensity of the peaks, 
the cj2 ^ correlations show a dip at the center ( Fig2c). 

As a final example we introduce quasiperiodicity by adding 
a potential Vj = 2A cos(27T7j + </>) where A is an amplitude 
and 7 is an irrational number which introduces competing pe- 
riodicities in the system. We choose 7 = (a/(5) — l)/2 and 
4> = 7r/4. The free-fermion Hamiltonian with this poten- 
tial exhibits localization-delocalization transition at A = 1 
Jl7ll . In our numerical study, 7 is replaced by a ratio of 
two Fibonacci numbers F n /F n+ i, = Fq = l,F n+ i = 
F n + F n ^x), which describe the best rational approximant, 
obtained by continued fraction expansion of 7 11711 . Here we 
treat the case of 7 = 55/89, L ~ 89 and N = 25. In the ex- 
tended phase (A < 1), the spectrum is effectively continuous 
and the eigenstates are of the Bloch type, while in the local- 
ized phase, the spectrum is point-like and the eigenstates are 
exponentially localized. Because in its ground state, all the 
lowest N single particle levels are occupied, the HCB system 
displays its metal-insulator transition in the density profile: in 
the extended phase, the on-site density varies smoothly be- 
tween sites, while in the localized phase it is discontinuous. 

The effects of quasiperiodic disorder on HCB are depicted 
in Fig. 3. We see in Fig. 3a that the localization transition de- 
stroys the cusp in the central peak of the momentum distribu- 
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FIG. 3: Quantum coherence as quasiperiodic disorder varies: A = 
O(black), 0.5(green), l(red), 2 (blue). Figures on the right show 
a blowup of the Fibonacci peaks which are barely visible on the 
left. Additionally, distinct harmonic peaks can be seen at 68, 76, 81. 
(Here q 2 = SI/ J = 0) 



tion, demonstrating the loss of first order quantum coherence. 
However, second order quantum coherence is preserved, as 
shown in the structure of the noise correlations, Fig. 3b. This 
is analogous to the changes in coherence associated with the 
Mott transition in the absence of disorder. Additional peaks 
reflecting quasiperiodic order of the system appear at the re- 
ciprocal lattice vectors =F-F n (Figs.3c,d). The localization 
transition destroys these peaks in the momentum distribution 
(3c), yet they remain in the noise correlations (3d). In other 
words, in the localized phase, no trace of the quasiperiodic 
order remains in first order correlations, but it is still seen 
clearly in second order correlations. The intensities of these 
second order Bragg peaks increase with disorder in the ex- 
tended phase and reach a maximum value at the onset to local- 
ization, after which they decrease. Furthermore, at the critical 
point as well as in the localized phase, we see additional struc- 
ture at various harmonics of two frequencies that underly the 
quasiperiodicity of the system. The harmonic peaks in the lo- 
calized phase reflect self-similar fluctuations of the wave func- 
tion 11811 . We hope that our studies of quasiperiodic quantum 
coherence (whose details will be published elsewhere) will 
stimulate experimental studies of such systems, for example 
in two-color optical superlattices fiUl . 



( G aa +1 G a b—l G a , 



M = 



b+1 



Gb-la+1 Gb-lb-1 Gb-lb+1 
Gb+la+1 Gb+lb-1 Gb+lb+1 



G a c \ 

Gb-ic 
Gb+u 



, x = 



/ Gba 

G a +la 

Gb-la 
G c +la 



V G c -ia+i G c _ib_i G c _ih+i .. G c _ic / 



Gbb-i Gbc+i ■■ Gbd-i Gbc \ 

Ga+lb-l Ga+lc+l ■■ G a +ld-l G a +lc 



Gb-lb-1 Gb-lc+1 
Gc+lb-1 G c +lc+l 



Gb-ld-1 Gb-lc 
Gc+ld-1 G c +lc 



\ Gda ■■ Gdb-i Gdc+i ■■ Gdd-i Gdc ) 



s = 



/ G aa G a b+1 ■■ G ac \ 
Gba Gbb+1 ■■ Gbc 



\ G c -ia Gc-lb+1 ■■ G c —lc ) 



Y = 



f Ga+la G a +lb Ga+la+1 ■■ G a +\b-l G a+lc+1 ■■ G a +ld-l \ 



Gb-la Gb-lb Gb-la+1 
Gc+la G c +lb Gc+la+1 



Gd-la Gd-lb Gd-la+1 
G C a G c b Gca+1 
\ Gda Gdb Gda+1 



Gb-lb-1 Gb-lc+1 
Gc+lb-1 Gc+lc+1 



Gd-lb-1 Gd-lc+1 
Gcb-1 G cc +1 
Gdb-1 Gdc+1 



Gb-ld-1 
Gc+ld-1 

Gd-id-i 

Gcd-1 

Gdd-i ) 



In the above matrices, dots indicate continuous variation of the indices. In the absence of dots, the indices are explicitly written 
and they may not change continuously. 
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